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Abstract 





This paper investigates a single server batch service queueing system with second optional service and 
reneging. The transient state probabilities of the queueing model are computed by using fourth order 
Runge Kutta method. Cost analysis is performed to determine the optimal values of the service rates of 
First Essential Service (FES) and Second Optional Service (SOS) simultaneously at the minimum total 
expected cost per unit time. Some important performance measures and numerical illustrations are 
provided in order to show the managerial intuitions of the model. 
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1.Introduction 


In queueing theory, a service process is among the 
important characteristics of a queueing system. 
Recently, most of the service systems adopted the 
provision of more than one service. For instance, 
in transport industry passengers get transport 
service from buses/airline and other services like 
food, soft drinks, etc. However, majority of them 
may prefer only main service while few of them 
may demand an optional service. Complete details 
of these situations have been presented by [1-5] 
and considered an M/G/1 queue with second 
optional service which assumed to _ follow 
exponential distribution while first essential 
service is a general distribution. [6-8] considered 
both SOS and FES having a general service time 
distribution. 

Balking and reneging are two impatient acts where 
a customer tends to leave without receiving the 
service due to longer waiting. Queuing problems 
with balking and reneging was studied by [2]. A 
single server queueing model with impatience has 
been presented by[9,10], where the customers lose 
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patience if the wait is more than the threshold they 
fixed. Later on, [11,12] studied a finite buffer 
multiple working vacation queues with balking, 
reneging and vacation interruption under N-Policy 
and obtained the solution for the steady state 
probabilities using recursive. 

The study of the transient solution of the 
performance measures of the queue system has 
been analysed by [14-18] etc. A time dependent 
solution of a single server queueing system with a 
Poisson input process has been considered by 
[19,20]. The transient behaviour of the infinite 
capacity M/G/1 modelwith batch arrivals and 
server vacations has been discussed by [21]. The 
transient behaviour analysis of an M/M/1/ 
N queue with working breakdowns and _ server 
vacations has been studied by [22]. 

From the above literature study, most of the 
heterogeneous server queueing systems are 
analysed in steady state. However, the steady state 
measures do not explain the real picture of the 
system behaviour, since they ignore the transient 
and the initial behaviour of a system [18]. In most 
of the practical applications, the state of the system 
experiences a change and such changes can be 
comprehended by the transient behaviour. 
Therefore, this article aims to study amM/M/1/1/ 
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Nqueue with SOS and reneging in transient state. 
The transient probabilities of the queueing model 
are obtained using fourth order Runge-Kutta 
method. 


2. Model Description 


We consider an M /M I/1/N queueing model 
with SOS and reneging customers in which the 
customers arrive according to a Poisson process 
with mean rate 1. The service times distribution of 
both FES and SOS are exponential with mean rates 
Ul, and Uz, respectively and the services are given 
in batches of size not more than b such that if the 
server finds the customers less or equal to b in the 
waiting queue, the server takes all of them in the 
batch for service, but if the server finds the 
customers more than b waiting in the queue, then 
the he takes a batch of size b while others remain 
waiting in the queue. The FES is required by all 
arriving customers and after completing FES, they 
may opt SOS with probability r or may leave the 
system with probability (1 — 1). The capacity of 
the system is finite (NV). A customer upon joining 
the queueing system and waiting in queue for some 
time, may get impatient and renege from the queue 
if his waiting time is beyond the expected time 
limit. Reneging times follow an _ exponential 
distribution with rate a. 


2.1 Formulation of Mathematical Model 


The M/M!*!1/1/N queue with SOS and reneging 
can be modeled by a two-dimensional continuous 
time Markov process {(N(t),/(t));t = 0}, where, 
N(t) is the number of customers in the queue, and 
J (t) is the server state with 


ne e if the server is providing FES 
iG) 2; if the server is providing SOS 


The state space of the Markov process is given by 
0 = {(n,i);n = 0;i = 1,2} 
The transient probabilities are defined as 


Pri(t) = PriN@) =n J(t) =i,0sn<N, 
i=1,2 


Here, P,,;(t) is the transient probability that there 
are n customers in the queue at timet and the 
server is providing FES (SOS) service, and Q(t) is 
the probability that the queue is empty at time t and 
the server is idle. Using Markov theory, 
thedifferential-difference equations of the model 
are as follows: 
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Q'(t) = -AQ(t) + 1 — eet aon 
+ [2 Po 2(t) ..-. pagdecaae wie, 


Post) = —A+ Wy) Poi) + 1 - 
T) iy Deed Pi (t) + be pan Pin (t) + a@Py4(t) + 


Pia glt) =-At w+ a)Pri(t) 
+ (1-1) py Patni (t) 
+ U2 Prip2(t) + @Prii1(t) 


+AP, ast), Lan 

SN Des seckieectnetienié Seb cee sede O 
P'na(t) = —A4F by + @) Pai () + @Pr4i1() 

+ APy-11(t), 

N-b+1<n<N-1......... 4 


P'yi(t) = —(p, + a)Py1(t) + APy-1,1(t) nee ees 5 


P'o2(t) = —A + 2) Po2(t) + rH1Po1(t) 
PP ya (0) sex ste ses sc cen etsy con cee'd sO 


P’y2(t) = —A+ pe + @)Py2(t) + APy_11(6) 
+ TpyPni(t) + nau 
lena N= 1, siatiperitaual 


P'y2(t) = —(U2 + @) Py 2(t) + 1Py_44(8) a5 
rU,Py,i(t), Mas Whe eeseesaesaeosisesddanesnassen ane smneess 8 


The above is an implicit set of differential 
equations whose analytical solution is a difficult 
task. However, we shall use the numerical 
technique to derive the approximate solutions 
assuming some initial conditions. 


3. Transient Solution of the Model 


In this section, we apply the fourth order Runge- 
Kutta method (R-K-4) to obtain the numerical 
solutions to the above equations. We assume that 
the system is empty at the starting point when the 
server is idle and the probability of the state of the 
initial system at time t = 0 are given as follows: 


Q(0) = 1,P,;(0) =0,i=1,2, O<n<N 


The set of differential difference equations (1) to 
(8) can be expressed as: 


xX 
crs f(t,x); x(t0) = x0 


The formula used for Runge-Kutta method of 
fourth order is given by: 
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k, + 2k. +2k3+k 
1 2 3 4) At 
6 
where; Atis the step size,andO < j < N -1 


ky = f (Gj, x; )At 


xy t1=x+| 


k 
= At 


At 
ky = Gt Soy + ) 


At ky 
kz = f(t a ae +> )At 


In our model, we consider x as P, ;(t)and f(t, x) 
aPyi(t) 


as 
dt 


4. Performance Measures 


In this section, we present the performance 
measures for anticipating the system behaviour. 
The various performance measures are as given 
below. 


(a) Expected number of customers in the queue 
during FES, SOS and the overall queue length 
are given, respectively by 


N 
LqFES(t) = ». nPy,1 (t), 
n=0 
N 
LgSOS(t) = ». nPp.o (t) 
n=0 


Lg(t) = LgFES(t) + LgSOS(t) 


(b) The average reneging rate is given by 
N N 


R.R(t) = ». aP, 1 (t) + > aP,, > (t) 
n=0 n=0 
(c) The blocking probability and the effective 
arrival rate are, respectively given by: 


Pblock(t) = Py y(t) + Py2(t), 
A(t) = A(1 — Pblock(t)) 


N-1 
=) aco + > (Pna(t) + Pn,2 ) 


(d) As a result of reneging, the Average Load 
(AL) on the server is given by 


AL(t) = A'(t) —R.R(t 
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Using Little’s Law, we can determine W(t) 
as follows: 


_ La 


Wal) = AL(t) 


q(t) 





4.1 Cost Analysis 


In this section, we develop the total expected cost 
function per unit time based on the steady state 
system performance measure when t > oo. Our 
objective is to determine the optimal value of 
service rates of FES and SOS ie., (u, and wy), 
respectively, so that the total cost function is 
minimized. Let us define the following cost 
elements: 


C, = Cost per unit time per customer present in 
the queue, 


C, = Cost per unit time when the server is serving 
in FES,C, = Cost per unit time when the server is 
serving in SOS, C, = Cost per unit time per lost 
customer due to reneging. 


From the definitions of each cost element listed 
above, the steady state expected cost function 
f (t41, H2) per customer per unit time can be written 
mathematically as: 


f(r be) = Cglg + Cy + Cov, + C,E[R]. 


The above expected cost function is_ highly 
complex and non-linear in terms of f (U4, U2)and 
the derivatives of the cost function are not easily 
available. Therefore, we solve the above 
optimization problem using Quadratic Fit Search 
Method (QFSM), since this method is highly 
suitable for the non-differentiable objective 
functions. 


Given the point € = (é,,€,), through QFSM, we 
fit a quadratic function that corresponds to 
functional values having a unique minimum, €7 = 
(é7 ,€7), for the given objective function f(€) = 
f (&1,&2). Quadratic fit makes use of this 
approximation to improvise the current 3-point 
pattern by replacing one of its pointswith 
approximate optimum é?%. The unique optimum 
€4of the quadratic function agreeing with f(€) at 
3 —point operation (&',€™, €" ) occurs at 
q 

2 fuk seis sl +7 EE os" 
~ APE = Pie EY = ol + FR — 6 
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Wheres! — (e)7 i,m — (e")" jst = (ey 
For a detailed algorithm, one may refer Rardin 
(1997). 

5. Numerical Investigation 


In this section we demonstrate the applicability of 
the Runge-Kutta method in transient state models 
using MATLAB - software. The numerical 
computations are presented in the form of tables 
and graphs. The values of the system parameters 


Table 1: Probabilities using Runge Kutta method 
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are assumed to beA = 2,u, = 3,u, = 2.5,a = 
0.1,r = 0.36,N = 10,b = 5. 


As observed in Table 1, the transient probabilities 
obtained by Runge-Kutta method tend to the 
corresponding steady state probabilities for large 
values of t. This confirms correctness, accuracy 
and efficiency of the model. 


Table 3: Optimum cost for different b and 
A(N = 10,r = 0.3) at steady state for t — 00 






































P,,,(t) Transient solution 
t=3 t= 10 t > © 
Q(t) 0.3471 0.3457 0.3457 
Po i(t) 0.2733 0.2731 0.2731 
Po 2(t) 0.0668 0.0668 0.0668 
Pui) 0.1096 0.1097 0.1097 
Py 2(t) 0.0555 0.0555 0.0556 
Pio1(t) 0.0000 0.0000 0.0000 
Pio2(t) 0.0001 0.0000 0.0000 
Total 1.0000 1.0000 1.0000 














Table 2: Variation in different measures of effectiveness with the change in arrival rate (A) 


p= lba=(. ar =O 


(5.54219, 3.15511, 38.9977) 0.607336 
(6.35025, 3.05648, 45.0453) (1.940152 

(7.08885 ,9.51125, Lie 1) 1.02275 
(8.6376,4.14159, 57.9119) 1.09593 


(0857459 
().LL0784 
0.115751 
0.120875 


().000493297 
().00212125 
(1.002 79075 
().00327201 


L327 ( 
2.38091 ( 
2.STOSS | 
3.00707 


17434 
94374 
Do0aL 


20426 


I 
I 
), 


3 
3 
4 
().3 


b= Ja = 0.5 


(2.83679,1.80106, 27.4659) 0.940152 
3 63091,2.0872, 31.3237) 0.782394 
(4.08024,2.1111, 34.8848) — 0.988905 
(4.50026,2.61151, 38.5909) 1.00624 
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0.181567 
(). 200987 
().211201 
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().0212125 394974 
().001 11562 0.337873 
(1.026094 2 0.300082 
(100204425 (306789 


ee) 
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Table 3 provides the sensitivity analysis for 
optimum servicep;,5, optimum cost f* and some 
performance measures. We observe that for fixed 
batch size the optimum service rates, 


cost, 
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arrival rate is increasing, which is in accordance 


with the practical situations. The trend is reversed 


Lq,R.R,Pblock, AL and W, are increasing as the 


Table 3 provides the sensitivity analysis for optimum servicey, uw, 


when we fix A and increase batch size value b. 






























































15 
Time 
Fig.1: Transient state probability for A = 2 
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a Lq Pblock X' Rr AL Wq 
1 0.1822 0.0000 0.9999 0.0413 0.9587 0.1900 
2 0.5937 0.0003 1.9994 0.0654 1.9340 0.3070 
3 1.1182 0.0029 2.9914 0.0783 2.9130 0.3839 
4 1.7160 0.0114 3.9543 0.0850 3.8693 0.4435 
5 nd INT 0.0284 4.8578 0.0879 4.7699 0.4939 
6 3.0034 0.0541 5.6754 0.0882 0871 0.5376 
7 3.6295 0.0868 6.3925 0.0869 6.3056 0.5756 
8 4.2135 0.1241 7.0070 0.0844 6.9226 0.6087 
9 4.7447 0.1638 7.5255 0.0813 74442 0.6374 
10 5.2199 0.2041 7.9590 0.0778 7.8812 0.6623 
” 
® 
CS 
oO 
wo 
2 
2 
ou 


Oo 
Oo 
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Probabilities 





Fig.2: Transient state probability for A = 5 
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Fig.3: Expected queue size versus time 
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Fig.4: Expected queue size versus time 
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Fig.5: Effect of a on the expected queue size with respect to time 
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Fig.6: Effect of r on the expected queue size with respect to time 


Figures land 2, plots the transient — state this is due to the fact that more customers demand 
probabilities versus time for different A = SOS. 
2.0and 5.0, respectively with a = 0.35. This 


graph is to show the sharp decrease in the idleness 
probability and the increase in FES probabilities as 
is the case when the rate of arrivals increases. 
Figure 3 shows the effect of arrival rate 2 on the 
expected queue size with respect to time. We 
observe that the expected number of customers in 
the queue gradually increases as A increases. In 
figure 4 one can observe that the expected queue 
size increases progressively from zero initially up 
to a certain value where it attains the steady state. 
In addition, we observed that the expected queue 
size in SOS is lower than in FES. Figures 5 and 6 
demonstrates the impact of a and r, respectively on 
the expected number of customers in the queue 
versus time. We can see that as @ increase the 
expected queue size decrease. Furthermore, we 
observed that initially the customers were patience 
with a variation of reneging rate, but as time goes 
on, some customers’ loss patience and decide to 
leave the queue. In figure 6 we observe that asr 
increases with time the queue size in both case 
SOS and FES increase. In addition, when r 
increase the expected queue size in SOS increase, 
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6. Practical Application of the model 


The model has numerous applications such as in 
hospital management systems, bank services, 
computer and communication networks, production 
system, etc. For example, in the hospital 
management systems the model can be applied to 
the situation where the outpatients request the 
appointment in the clinic system. The clinic officer 
monitors the length of booking windows for 
appointment of the outpatients, but since in most 
cases there iS limited number of 
physicians/doctors,it leads to the unbalanced ratio 
between the number of outpatients and that of the 
physicians/doctors. This situation leads to an 
increase in the length of the booking window and 
brings the long waiting for appointments. 
Therefore, in order to shorten outpatients waiting 
time, we limit the length of the booking window 
and we assume the system has a limited slots 
capacity N during a limited length of booking 
window. The system capacity N can be divided 


: . N 
into x equal amount of slots bi.e.(b = . ). The 
slots b are termed as a maximum size of a batch of 
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the outpatients appointment and x is equal to the 
number of sessions/periods. In practice, the new 
requests for appointments come individually and 
the service can be provided from single patients up 
to the maximum batch size b of the patients. And 
after the first service the patients may opt for 
another service provided by a clinic. On the other 
hand, some patients cancel the appointment and 
those patients are termed as reneging patients. 


Conclusion 


In this paper, we present M/M!!/1/N queue with 
second optional service and reneging in transient 
state. The fourth order Runge Kutta method is used 
to determine the transient solution of the queueing 
model and obtained the transient probabilities of 
the model. We developed the cost analysis to 
determine the optimal service rates of FES and 
SOS. Furthermore, we _ presented some 
performance measures along with the numerical 
results to show the effect of various system 
parameters on the transient behaviour of the model. 


References 


[1] Al-Jararah, J andMadan, K.C. (2003). An 
M/G/1 queue with second optional service 
with general service time distribution, 
International Journal of Information and 
Management Sciences, 14, 47—56. 


[2] Ancker, C.J andGafarian A.V. (1963). 
Some queueing problems with balking and 
reneging, Operations Research, 11, 88-100. 


[3] Ke. J.C..Wu, C.H and Wen Lea Pearn, 
W.L. (2013). Capacity multi server queue 
with second optional service channel, 
International Conference on Mathematical 
Modelling in Physical Sciences, Conference 
Series. 


[4] Madhan, K.C. (2000). An M/G/1 queue 
with second optional service, Queueing 
Systems, 34,3746. 


[5] Madhan, K.C.(1994). An M/G/1 queueing 
system with additional optional service and 
no waiting capacity, Microelectronics and 
Reliability, 34,521—527. 


[6] Manoharan, P and Sankarasasi, K. (2015). 
An M/G/1 reneging queueing system with 
second optional service and second optional 
vacation, Applied Mathematical Sciences, 
9, 3313-3325 


International Research Journal on Advanced Science Hub (IRJASH) 


Volume 02 Issue 10 October 2020 


[7] Manoharan, Pand Sankarasasi, K.(2015). 
An M/G/1 feedback queueing system with 
second optional service and second optional 


vacation, International Journal of Fuzzy 
Mathematical Archieve, 8,101— 113. 


[8] Medhi, J. (2002). A single server Poisson 
input queue with second optional channel, 
Queueing Systems, 42, 239-242. 


[9] Wang, J. (2004). An M/G/1 queue with 
second optional service and _ server 
breakdowns, Computers and Mathematics 
with Applications, 47, 1713-1723. 


[10]Choudhury, A. (2008). Impatience in 
single server queuing model, American 
Journal of Mathematical and Management 
Sciences, 28, 177-211. 


[11]Neuts, M.F. (1967). A general class of bulk 
queues with poisson input, The Annals of 
Mathematical Statistics, 38, 759-770. 


[12]Laxmi, P. V, Goswami, V and Jyothsna, K. 
(2013). Optimization of balking and 
reneging queue with vacation interruption 
under N-Policy, Journal of Optimization, 
Vol. 2013. 


[13]Nagarajan, R and Kurose, J. F. (1992). On 
defining, computing and _ guaranteeing 
quality of service in high speed networks, 
Proceeding of INFOCOM92, Florence, 
Italy, 2013,2016—2025. 


[14]Holt, B. V and Blondia, C. (2005). 
Approximated transient queue length and 
waiting time distribution via steady state 
analysis, Stochastic models, 21, 725-744. 


[15]Hofkens, T, Spacy, K and Blondia, C. 
(2004). Transient Analysis of _ the 
DBMAP/G/1 queue with an application to 
the dimensioning of video playout buffer 
for VBR __ Traffic, Proceeding of 
Networking, Athens Greece. 


[16]Lucantoni, D. M, Choudhury, G. L, and 
Witt, W. (1994). The transient BMAP/P 
H/1 queue, Stochastic Models 10 461-478. 


[17]Bohm, W andMohanty, S. G.(1994). 
Transient analysis of queues’ with 
heterogeneous arrivals, Queueing Systems, 
18, 27-45. 


[18]Whitt, W.(1983). Untold horrors of the 
waiting room, what the equilibrium 


www.rspsciencehub.com 


distribution will never tell about the queue 
length process, Management Science, 29, 
395408. 


[19]Parathasavathy, P. P. (1987). A transient 
solution of an M/M/1 queue: A_ simple 
approach, Advanced Applied Probability, 
19, 997-998. 


[20]Bai, Y, Song, J  andZhang, G. 
(2015).Design appointment in outpatient 
department with bulk service, TEEE 
International Conference on Automation 
Science and Engineering (CASE), 746-751. 


[21]Jacob, M. J. and Mathusoodanan, T. 
P.(1988). Transient solution of an M/G/1 
queueing system with group arrivals and 
vacations to the server, OPSEARCH, 25, 
279-284. 


[22]Yang, D. Y andYing, Y. W. (2014). 
Transient behaviour analysis of a finite 
capacity queue with working breakdowns 
and server vacations, Proceedings of the 
International Multi Conference of 
Engineers and Computers Scientists, Hong 
Kong,2, 12-14. 


International Research Journal on Advanced Science Hub (IRJASH) 


Volume 02 Issue 10 October 2020 


